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Abstract 



The possibility of trapped modes of gravitational waves appearing in stars with R > 3M is considered. 
It is shown that the restriction to R < 3M in previous studies of trapped modes, using uniform density 
models, is not essential. Scattering potentials are computed for another family of analytic stellar models 
showing the appearance of a deep potential well for one model with R > 3M. However, the provided 
' example, although having a more realistic equation of state in the sense that faound < oo, is unstable. 

£C) | On the other hand it is also shown that for some stable models belonging to the same family but having 

R < 3M, the well is significantly deeper than that of the uniform density stars. Whether there are 
0^ ' physically realistic equations of state which allow stable configurations with trapped modes therefore 

remains an open problem. 

00 
On 

O ! 1 Introduction 

| | Examples of trapped modes of gravitational waves in compact stars were first given by Chandrasekhar and 
(DJO' Ferrari [0 and calculations were subsequently also carried out by other authors @, |l f|] . The fundamental 
reason behind the occurrence of the trapped gravity wave modes is the stretching of the geometry by the 
strong gravitational field leading to a bell-like geometrical structure inside the star. This phenomenon is 
most clearly illustrated using the concept of the optical geometry as developed by Abramowicz and coworkers 
(see H and references therein). The optical geometry of the vacuum Schwarzschild metric develops a neck 
precisely at R = 3M implying that for stars with R < 3M there will be a family of closed null geodesies in 
the stellar interior. It is natural to associate this behavior with the trapping of certain modes of gravitational 
radiation although the relation between the trapping and appearance of the neck in the optical geometry is 
only approximate. The optical geometry is useful not only for pedagogical purposes but can also be used 
to motivate an estimate of the eigenfrequencies of the resonances J5| . In previous studies it has often been 
assumed that the appearance of a neck (and the consequent trapping of gravity waves) is only possible if the 
star is ultracompact, that is the compactness, (3 := M/R, must lie in the range | < (5 < | where the upper 
bound | is Buchdahl's limit Q representing the maximum compactness for any static star for which the 
energy density is decrasing outwards. The compactness is usually given in terms of the inverse compactness 
a := R/M which we will refer to as the tenuity. The trapped modes found in [jjj occur for tenuities in the 
range 2.25 < a < 3. Realistic neutron stars are believed to have tenuities in the range 3 < a < 11 so they 
are at most marginally ultracompact in this sense Q. However, as will be shown in this letter, trapped 
modes may occur in stars with a > 3. This opens up the possibility for real neutron stars to exhibit gravity 
wave trapping. In view of this result it seems like a good idea to reserve the notion of ultracompactness 
for stars which have a neck in their optical geometry and consequently a family of closed null geodesies in 
their interior. Ultracompact stars would then be expected to exhibit gravity wave resonances as well. As 
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will become more clear later ultracompactness in this sense really applies to the stellar core rather than the 
entire star. Although the definition of ultracompactness given here is unambiguous it is more difficult to 
calculate in practice. In the concluding remarks we will touch upon possible rules of thumb criteria which 
could be used as a rough estimate of compactness. 

It is not difficult to understand why stars with R > 3M could have an ultracompact core. The key is the 
behavior of the equation of state at low pressure. Consider a uniform density model with radius less than 
3M. Now replace a thin shell (its mass should be finite but be only a small fraction of the total mass) at the 
surface with some material with a soft equation state (for example a polytrope) such that the total mass of 
the star remains the same. In physical terms we can think of this process as giving the star an atmosphere 
by transforming some of the matter near its surface. Clearly the gravitational field in the core is the same 
as it was before. However, the radius will depend sensitively on the equation of state of the atmosphere. In 
fact it can be made arbitrarily large for example by letting the atmosphere be a polytrope of index 5 — e 
where e <C 1. Another alternative would be to replace the slice by an envelope which, like the core matter, is 
of uniform density but satisfying /^envelope < Pcorc- Such double layer uniform density models were recently 
considered by Lindblom || to discuss phase transitions in compact stellar models. The radius could then be 
made aribtrarily large by letting the quotient ^envelope/ Pcoro be sufficiently small. 

Although the argument given above should be sufficient to establish the existence of trapped gravity 
wave modes for stellar models with R > 3M, there remains some critical issues concerning the realization 
of such models in nature. One such issue is the question of causality. Of course, already the unform density 
models are unrealistic in this sense having an infinite speed of sound. A second issue is that of stability. The 
absence of a local mass maximum in the uniform density models shows that they are in fact stable. In this 
letter we will use the generalized Buchdahl n = 5 polytrope (GB5) family of exact models |f^| to illustrate 
the new possibilities which occur when one considers softer equations of state. This family generalizes the 
original Buchdahl solution which behaves as a polytrope of index 5 at low pressure. The generalized models, 
however, have an equation of state which is liquid-like at low pressure in the sense of having p s > ("s" 
denoting the value at the stellar surface). 

2 Stellar models 

The metric of static spherically symmetric (SSS) models is usually given in the Schwarzschild form 

g = -e 2v dt 2 + e 2A dr 2 + r 2 (d(9 2 + sin 2 6»d0 2 ) . (1) 

For our purposes we also need to write the metric of a SSS system in a general radial gauge as 

g = -Y 2 dt 2 + N 2 dx 2 + S 2 {d9 2 + sin 2 (9d0 2 ) , (2) 

where Y, N and S are functions of the radial variable x. The Schwarzschild radial variable is then given 
by the relation r = S(x). Before proceeding we need to deal with a possible source of confusion relating 
to the metrics ([!]) and (g). The time coordinate is a priori only defined up to a scaling and a translation. 
The scaling gauge can be fixed by the requirement that the time coordinate should correspond to the proper 
time of a static observer at infinity. We shall refer to this gauge as the proper time gauge. This gauge is 
usually but not always imposed when writing down the metric of exact solutions. It is assumed here that 
the metrics (Q) and (|^) refer to the proper time gauge. Correspondingly the formulas given below are also 
given in this gauge. However, since exact solutions are not automatically given in the proper time gauge it 
is useful to write down the relevant transformation formula for a metric written in a general time gauge. To 
do that we first note that for the Schwarzshild exterior metric (as usual expressed in the proper time gauge) 
k := yl — 2M/ R = e" s = Y s , where the subscript s denotes the surface of the star. Therefore we must have 
e"" = k for the stellar model. Now let t be an arbitrary time coordinate and v (or Y) the corresponding 
metric functions. Then the required relations are 

e v = ke p -^ , Y = kY^Y , di=kY~ 1 dt . (3) 
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The gravity wave modes discussed in |lj are equivalent to non-radial axial (i.e. odd parity) perturbation 
modes of SSS fields. Such axial perturbations do not couple to fluid motions in the star a fact which accounts 
for their alternative interpretation as gravitational wave modes. The axial modes with frequency u> and mode 
number I > 2 are governed by the equation Jll| 



d 2 Z 
dr? 



VZ = u> 2 Z , 



(4) 



where the potential V = Vi + Vd is here written formally as a sum of a centrifugal and a dynamical part (cf. 
in the form 



V, =l(l + l)e 2v r- 2 , 

V d = e 2tJ [r- 1 (A, r - u, r + 2^) e~ 2X - 2r~ 2 } , 
and r* is the tortoise radial variable defined by 

dr» = e x - u dr = Y~ x Ndx . 

To express the potential in a general radial gauge we use the relations e v — Y , e A = NS'^ 1 and 

X, r = N- 1 N , S'~ 1 -S'~ 2 S" , 

where the primes denote differentiation with respect to x. The potentials then become 
Vi = 1(1 + 1)Y 2 S- 2 , 



(5) 



(6) 



(7) 



Ki = Y 2 N- Z S 



1 rr-2 



(N^N' - Y^Y') SS' - SS" + 2S' 2 - 2N 2 



Using the Einstein equations the dynamical part of the potential can be written in the form 

6m(r) 



V d = Ge l 



4n(p-p) - 



where 



m(r) = 47r / pr 2 dr = <iir pS 2 S'dx 
Jo Jx c 



(8) 



(9) 



(10) 



is the mass within radius r. We are using units in which c = 1 but keep the gravitational constant, G, for 
convenience in some formulas. Geometric units can be obtained by setting G = 1. 

The exterior Schwarzschild solution 

In this case 



,2i/ 



2GM 



(11) 



and m(r) = M leading to 



Vi = 



1(1 + 1)(1 - 2GM/r) 



Vd = - 



6GM(1 - 2GM/r) 



(12) 



V 



(1 - 2GMr)[l(l + 1) - 6GM/V] 
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The interior Schwarzschild solution 



Schwarzschild's uniform density model is characterized by 



e"=§(3fc-i0, e x = y~ l 



, 2G/3 3 r 2 3/? 3 , — 



(13) 



where /3 := M/R. The potentials then become 



^ + l)(3fc-;/) 2 3G/3 3 fc(3fc-y) 
= ^2 . V d = -47rGkpe = — 2 . (14) 

The GB5 interior solutions 

The family of exact interior solutions which we focus on this letter is the GB5 family || [To) given by (using 
a non-proper time gauge) 

f= T~rf ' N = S=(T + X) 2 , (15) 

where 



T(x) = J , A . b .. , JC(aj) = W— '— . (1(>) 

v ; ycosh(x-A) v ; V cosh a; 

The constants b and A characterize the equation of state while A is the single nontrivial integration constant 
appearing in all SSS models. The equation of state can be written in the form 

a(u 6 -A 6 ) 3a(u 5 + A 6 ) 

(17) 
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y {\ + uf(\-u) ' (1+u) 

where a := l/(87rG6 2 ), u := X/T and A := u s . Defining x '■= u c = A 3 / 2 e A / 2 ("c" denoting value at the 
center) we have 

0<A = u s <tt<it c = x<l- (18) 



The two parameters a and A which characterize the equation of state in (|17j) can be interpreted as a scaling 
and a stiffness parameter respectively. The scaling parameter a just represents a change of overall scale. All 
other physical characteristics in the model are unaffected by changes in a which can be any positive number. 
It is convenient to replace \ by another parameter /j, defined by the relation \ — A + M — A/i. In that way the 
A/u-section of the parameter space is exactly the unit square, 0<A<1, < /i < 1 (see figure 1). However, 
in order not to complicate the formulas unnecessarily we keep using \ but think of it as a function of A and 
p. For our purposes it is also useful to replace a by the mass of the star. Expressing a (or b) and A in terms 
of M, A and x we have 



where 

(X 2 + A 3 )(x 2 -A 3 ) 3 M 



4 x (x 2 - A 2 ) 3 / 2 (x 2 - A 4 ) 3 / 2 



(20) 
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In calculations we wish to use the set (A, a, M) as input parameters to specify the stellar model. The tenuity 
is given by the expression 

(1 + A) 2 ( X 2 -A 3 ) 2 (21) 



2A(x 2 -A 2 )(x 2 -A4) ■ 
In order to use a as an input parameter we solve this equation for \ which yields 

X - X {T—() ' C -(l + A>-(l + A)2 ' (22) 

In the radial gauge N — S the expression for Vd reduces to 

V d = y- 2 S- 4 (35' 2 - Y^Y'SS' - SS" - 2S 2 ) . (23) 

Inserting the GB5 functions Y and S in this expression gives the potential in an explicit but complicated 
form and we do not write it down here. 

3 Discussion 

We now consider the new possibilities which occur when using an equation of state with non-uniform density 
using the GB5 family as a theoretical laboratory. The parameter space of the GB5 models is shown in 
figure 1. 

[Figure 1 about here.] 

One of the two models marked in figure 1 has a = 5. The scattering potential for that model is shown in 
figure 2. 

[Figure 2 about here.] 

This clearly illustrates the fact that the potential may have a minimum in the stellar interior even though 
the outer parts of the star extend to regions well beyond a = 3. We also mention without proof that this 
model admits a family of closed null geodesies in its interior. It may be objected that the model is unstable 
(as indicated in figure 1) and that this result therefore has little physical relevance. However, the instability 
is closely connected with the softness of the equation of state. Taking instead the double layer uniform 
density models mentioned in the introduction it should be possible to provide examples of stable models 
having a potential with a minimum in the interior. A second comment we wish to make on this issue is 
that resonances in unstable models may in principle be important in gravitational collapse situations where 
short-lived unstable equilibrium states could perhaps form en route to the final collapse. 

The second model indicated in figure 1 has a = 2.4 and lies in the stable region of the parameter space. 
The corresponding potential is plotted in figure 3. 

[Figure 3 about here.] 

The phenomenon we wish to illustrate here is that the GB5 potential has a significantly deeper potential 
well than a uniform density model with the same mass and R/M ratio. The quasi-normal modes of the 
uniform density model with a = 2.4 were calculated in M. It would be interesting to calculate the modes 
for the GB5 model. The deeper minimum is an indication of longer damping times compared to the uniform 
density case. 

The question of whether realistic stellar models can be ultracompact (in the sense defined in this letter) 
remains open. In M], Iyer, Vishveshwara and Dhurandar searched for stable and causal models satisfying 
R < 3M. In view of the results given in the present work it would be more relevant to look for stable and 
causal models which are ultracompact in the sense of having a family of closed null geodesies in the stellar 
interior. It would be very useful to have a simple rough criterion of ultracompactness expressed in terms of 
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a dimensionless combination of easily computable quantities. Examples of criteria of compactness include 
the central redshift (defined as the redshift of a hypothetical speed of light signal sent from the center of 
the star and received by a static observer at infinity) and the central 4-dimensional curvature, for example 
(R a ^ lS R a 0-ys)c or (R a P R a p) c . The redshift is already dimensionless while the curvature measures need to 
be properly normalized for example by multiplying by a power of the total mass. However, it is not clear 
whether any of these measures, either by themselves or by taking combinations, could serve as criteria for 
ultracompactness. 
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List of Figures 



1 The A/i parameter space of the GB5 family. The pressure and energy density satisfy the physical 
requirements p > 0, p > and dp /dp > 0. Contours are shown for some values of the tenuity a — 
R/M. Specifying the value of A fixes the equation of state up to a scaling. The corresponding sequence 
of stellar models can then be followed along the A = constant line starting from the Newtonian limit at 
p, — 0. The mass increases along such a sequence up to the dashed curve after which it decreases. The 
dashed curve therefore represents the transition to unstable models. The stable region is consequently 

located to the lower right of the stability limit. The dotted curve is the causal limit. A model is 

classified as causal if the equation of state satisfies dp /dp < 1 throughout the interior of the star. Foi 

acausal models the speed of sound as defined by ti soun d i= \J dp/dp is therefore larger than the speed 

of light in some part of the stellar interior. The causal GB5 systems are located at the lower left side 

of the dotted curve in parameter space. The two marks in the figure are the points which correspond 
to the models for which potential plots are given in this letter 

2 The GB5 potential for the model marked on the a — 5 contour in figure 1 is plotted for / = 2 along 
with the exterior potential. For comparison the potential of the uniform density model with the same 
mass and radius is shown as the dashed curve 

3 The potential for the second GB5 model marked in figure 1 is plotted for 1 — 2. As in figure 2 the 
dashed curve is the potential of the corresponding uniform density model 
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Figure 1: The A/i parameter space of the GB5 family. The pressure and energy 
density satisfy the physical requirements p > 0, p > and dp/dp > 0. Contours are 
shown for some values of the tenuity a = R/M. Specifying the value of A fixes the 
equation of state up to a scaling. The corresponding sequence of stellar models can 
then be followed along the A = constant line starting from the Newtonian limit at 
a = 0. The mass increases along such a sequence up to the dashed curve after which 
it decreases. The dashed curve therefore represents the transition to unstable models. 
The stable region is consequently located to the lower right of the stability limit. The 
dotted curve is the causal limit. A model is classified as causal if the equation of 
state satisfies dp/dp < 1 throughout the interior of the star. For acausal models the 
speed of sound as defined by « SOU nd := \J dp/dp is therefore larger than the speed of 
light in some part of the stellar interior. The causal GB5 systems are located at the 
lower left side of the dotted curve in parameter space. The two marks in the figure 
are the points which correspond to the models for which potential plots are given in 
this letter. 
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R/M = 5, X = 0.1294 
0.3n < ; 
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Figure 2: The GB5 potential for the model marked on the a = 5 contour in figure 1 
is plotted for / = 2 along with the exterior potential. For comparison the potential 
of the uniform density model with the same mass and radius is shown as the dashed 
curve. 
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R/M = 2.4, X = 0.61 
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Figure 3: The potential for the second GB5 model marked in figure 1 is plotted for 
1 = 2. As in figure 2 the dashed curve is the potential of the corresponding uniform 
density model. 
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